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ABSTRACT 

The Very Large Array has been used at five frequencies to study the structure 
and linear polarization of SS433 on scales as small as ~0.1" ~ 500 AU. Each jet 
consists of a sharp, curving ridge- line at the leading edge, plus significant trailing off-jet 
emission, showing that they are enveloped by diffuse relativistic plasma. No kinematic 
model with constant jet speed fits our images on all scales, but they are consistent 
with variations in jet speed of around ±10% around the optical value. Our images 
show continuous jets with bright components occurring simultaneously in the two jets 
roughly every 35 days. When corrected for projection effects and Doppler boosting, the 
intensities of the two jets are intrinsically very similar. Fractional linear polarization 
up to 20% is present along the ridge-lines, while the core is essentially unpolarized. 
The rotation measures are consistent with a foreground screen with RM ~ +110 rad 
m -2 plus a larger, asymmetrical contribution close to the core. The apparent magnetic 
fields in the jets are roughly aligned with the ridge-lines in most but not all of each 
jet. The jet is more highly polarized between the components than in the components 
themselves, suggesting that the mechanism that creates them compresses a tangled part 
of the magnetic field into a partially-ordered transverse layer. The off-jet emission is 
remarkably highly polarized, with m ~ 50% in places, suggesting large-scale order of the 
magnetic field surrounding the jets. This polarized signal may confuse the determination 
of magnetic field orientations in the jets themselves. However, the images are consistent 
with a jet magnetic field that is everywhere parallel to the helices. 

Subject headings: binaries: close — magnetic fields — polarization — radio continuum: 
stars — stars: individual (SS 433) 



1. INTRODUCTION 



The galactic binary source SS 433, consisting of a stellar-mass black hole in cl ose orbit about an 
early-type star, is the prototypical "micro-quasar," a miniature version of an AGN (IMirabel &: Rodriguez 
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19991 ), Two relativistic jets emerge from opposite sides of the compact object at speeds of 0.26 c. 
Modeling of the optical spectr um shows that the jet system precesses wi t h a period of 162 days 
abou t a cone of half-angle 20° (jAbell fc Margonlll979l ; Fabian fc Reeslll979l ; lMilgromlll979l ; iMargon 
19841 ), Imaging by the Very Large Arr ay (VLA) confi r med this picture, as the radio images showed 
helical jets on both sides of the source ([Spencerlll979l ; iHiellming &: Johnstonlll981allbr). Higher res 



oluti o n ima ges made by VLBI reveal the structure down to a scale of a few AU (jVermeulen et al 
1987Lll99ah . 



Arcsecond-scale imaging with the VLA at 5 GHz showed that the radio jets are linearly po- 
larized at levels u p to tens of percent, and tha t the polarized structure varies with position and 
precession phase (jHiellming &: Johnston! Il981al ), Hjellming and Johnston inferred that the un- 
derlying magnetic fields were roughly aligned with ballistic (i.e., radial) motions in the jets, but 
acknowledged that in the absence of any detail ed information on Faraday rotation in SS433, this 

J198l|) also reported strong and variable polarization in 



conclusion was tentative. 



Gilmore et al 



the jets at two epochs. IStirling et al.1 (J2004J) observed SS 433 at multiple frequencies, and concluded 
that the jet magne tic fields lay along the spiral of the jet. At VLBI resolutions SS433 shows no 
linear polarization (|Paragill2002l ). There is some circular polar ization when observed at 1-9 GHz 
by the Australia Telescope Compact Array (jFender et al.ll2000l ). 



In this paper we use the VLA to study the structure, linear polarization, and Faraday rotation 
of SS433 on scales as small as ~ 0.1" ~ 500 AU. In $2] we describe the observations and data 
reduction. In $3]we determine the properties of the jets, find the underlying magnetic field orien- 
tations, and discuss the physical implications of these results. Our conclusions are summarized in 



2. OBSERVATIONS AND DATA REDUCTION 

Interferometer data for three ~4 hour periods over 1998 March 5-7 (MJD 10877-10879) were 
obtained from the VLA data archive. The array was in the A configuration, with 26 working 
antennas; the frequencies used were 1.415, 1.665, 4.800, and 14.94 GHz on March 5, and 1.415 , 1.465, 



4.800 , and 14.94 GHz over March 6-7. The data were edited and phase calibrated in AIPS (INRAO 



20041 ) using the nearby unresolved source J1950+081, and amplitude calibrated against 3C286. 
Instrumental polarization contamination was determined using the phase calibrator, and removed 
from the data in the standard manner, and electric vector position angles calibrated against 3C 286 



The d ata were then imaged in total intensity and linear polarization using DIFMAP ([Shepherd et al 



19941 ). utilizing several cycles of phase and amplitude self-calibration. Natural weighting was used 



at 14.94 GHz, and uniform weighting at the lower frequencies. Images in Stokes Q and U were 
combined to make images of polarized intensity p = \J Q 2 + U 2 , electric vector polarization angle 
X = (1/2) arctan (U/Q), and fractional linear polarization m = p/I. 



Images were made at four frequencies for each of the three days, and they were consistent 
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with each other, and with the small proper motions expected in the jets over this time interval 
(fi ~ 9 masd -1 ). The core varied by a few percent at every frequency at which we had more 
than one day's data. To optimize the signal-to-noise ratios in the jets, we aligned the images on 
the core, and simply averaged the multiple images at each of 1.415, 1.465, 4.800, and 14.94 GHz. 
Smearing of the jet features over the span of the observations amounts to only 15% of a beamwidth 
at 15 GHz, and does not result in false structures, either in total intensity or in polarization. Along 
with the single image at 1.665 GHz, the results are presented in Figures [THSJ We also created a 
14.94 GHz image convolved to the resolution of the 4.80 GHz images (not shown). In each case the 
top panel shows contours of total intensity, the middle panel contours of linearly polarized intensity 
with electric vector position angle ticks, and the lower panel fractional polarization in grey scale 
(color in the electronic edition) over every-other total intensity contour. The images at 4.8 and 
15 GHz are consistent with pub li shed images at different e pochs and precession phases made by 
others (e.g., IStirling et al.1 (|2004l ); iBlundell Bowlerl (|2004l )). In particular, the helical structure 
out to 3" from the core in Figure [H is real, and is not an artifact due to the sidelobe structure of 
the observations. This is supported by the agreement of the images and the kinematic model (see 
below) . 

On each image we have plotted model jets, with positions calculated from the kinematic model 
of Hjellming and Johnston (1981b) , in which mater i al is e mitted as identical oppositely-directed 
narrow jets. In the convention of IVermeulen et al.l (|1993l ) (phase is zero when the east jet has 
maximum redshift), the precession phase at this epoch is (f) = 0.89. Shown are predicted positions 
of material emitted at five-day intervals; material that is oncoming is marked with "+," while 
retrea ting material indicated with "o." We adopted the optical ephemeris of lEikenberry et al. 
(|200ll ). combined with a position angle of the jet preces sion axis of PA = 98.2° determined from a 



sequence of VLB A images made in the su mmer of 2003 (Mioduszewskill2006l ): see Table HJ Several 
distances have been pr oposed for SS433. Mioduszewski ( 20od ) found a dist ance of d = 5 . kpc , 
Vermeulen et al.l (|1993l ) concluded from VLBI imaging th at d = 4.85zfc0 . 2 kpcjDubner at al.1 (|1998l ) 
derived a distance of 3 kpc from H I absorption, while ISthiing et al.l (120021) found a distance o f 
4.8 kpc, with jet deceleration at the rate of —0.02c per precession period^ Blundell & Bowler (2004) 
conclu ded from a deep VLA image at 5 GHz that d = 5.5 kpc. Recently. lLockman. Blundell. Goss 
(|2007l ) hav e re-examined the H I dist ance, and conclude that it is in agreement with the kinematic 
distance of lBlundell Sz Bowlerl (120041 ) . In Figure [6] we plot our 4.80 and 14.94 GHz images against 
three kinematic models, (a) one with d = 5.0 kpc and v = 0.2647c, (b) one with d = 5.5 kpc and 
v = 0.2647c, and (c) one with d = 5.5 kpc and v = 0.29c. The fit to the inner parts of both jets at 
14.94 GHz seems to favor eithe r (a) or (c), but model (b ) fits the larger-scale 4.80 GHz image. I n 
accord with the conclusions of IBlundell &: Bowlerl (|2004l ) and lLockman. Blundell. Gosd (|2007l ) , 
we adopt d = 5.5 kpc, and interpret the discrepancy in fitting the 14.94 GHz images as a sign 
th at synchronized varia t ions in velocities of about ±10% tak e place in the two jet s, as suggested 
bv lstirling etaD JsQoJ, blundell Bowlerl (fcooi bood ). and blundell et al.l (hoOTh . Since we are 
interested primarily in the smaller scale s tructure, in all fi gures we use d = 5.5 kpc and v = 0.29c. 
We also ignore the 6.5 d nodding motion (jKatz et al.lll982l ) because it is small on the scales we are 
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considering. 



3. PROPERTIES OF THE SS 433 JETS 
3.1. Structure 

The total intensity structure of SS 433 at 14.94 GHz (Figure [5|) shows seven distinct bright 
spots (the core plus six "components") that we denote El, E2, E3, C, W3, W2, and Wl, as labeled 
on Figure El and as indicated by dots on the other images. We chose those locations to study the 
properties of the jets in SS 433. The 14.94 GHz total intensity image was fit to the components 
using the AIPS task JMFIT, and the resulting positions used at all frequencies. We measured /, 
p, and x a t each frequency at each component by reading image values at its location. The results 
for 4.80 and 14.94 GHz are listed in Table Since we do not resolve the components cleanly at 
frequencies other than 14.94 GHz, the corresponding values should be regarded as estimates. The 
polarization detected at C is sufficiently small that we do not consider it further. 

We can use the kinematic model to estimate the birth epochs of the six components. Assuming 
that the precession rate of the jets does not vary while the velocities of ejections may, for components 
E1-E3 and W1-W3 we used their position angles with respect to the core to find their ages. For 
E2, E3, W2, & W3 the model points for (3 = 0.29 and d = 5.5 pc lie on the observed components. 
For El and Wl, the model predicts that the components should be further from the core than 
they are. We attribute this to a smaller ejection speed of j3 ~ 0.26 (essentially the optical value) 
at the births of these components, which is also the speed needed to fit the larger-scale structure 
at 4.80 GHz at the locations of the first loops in the jet. The resulting ages are given in Table EJ 
In addition, Figure [5] contains hints of two components born later than E3 and W3. This is 
better seen in a 14.94 GHz image from which the core has been removed, shown in Figure [71 We 
denote these tentative components E4 and W4. These components are separated from the core 
by only 0.7-0.9 beamwidths and are much fainter than the core, so we do not include them in 
the statistical analysis; however, they are included in Figures IHHTOl Because their position angles 
are quite uncertain, the ages of E4 and W4 were estimated using their distances from the core, 
with the results r(E4) ~ 12 days and r(W4) ~ 11 days. The model ages of the east and west 
components of each of the four pairs are plotted against each other in Figure [HJ It is clear that 
within the errors, the east and west components of each pair were emitted simultaneously. The 
intervals between emission of the components on a given side range between 27 and 41 days, with 
a mean of 35 ± 3 days, but the events do not appear to be strictly periodic. Including components 



E4 an d W4 does not alter these conclusions. We note that the 15 GHz VLA image of lStirling et al 



([2004]), made at precession phase close to ours {<p = 0.97), contains jet peaks w ith similar intensitie s 



and spacings to the ones in Figure [5] (see their figure 8). It is interesting that iFiedler et al.1 (119871 ) 
found a similar characteristic timescale for flaring in SS433 of 25 ± 10 days. 



At 14.94 GHz, we are able to resolve across the jets, and see that the bright ridge-lines are 
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"leading edges" of both jets, and are significantly sharper than the trailing edges (Figure [5]). There 
is a significant off-jet radiation, most notably north of the jet ridgeline to the east, and south of 
the ridgeline to the west. Most, but not all, of this emission is within the cone of precession of the 
jets. This extended plasma may be shed from the jets, or ambient medium in which relativisitic 
electrons are energized by interaction with the jets. 

It is interesting to determine whether or not the east and west jets are intrinsically identical. 
We can predict the run of brightness down the jets using the kinematic model, including both 
projection and Doppler effects. In particular, we can calculate the expected ratio of east to west jet 
brightness, and compare that to the ratios observed in the pairs of peaks El &: Wl, E2 &: W2, and 
E3 & W3. We resolve the components from each other and from the core only at 14.94 GHz, so we 
consider only that frequency. In all cases we adopt d = 5.5 kpc and a velocity of v = 0.29c. First, 
assuming that it is optically thin, the observed brightness of the jet depends on projection effects, 
essentially how much of it is included in a resolution element. We predicted the local density of jet 
material on the sky from the kinematic model by simply counting the number of model points per 
beam area along the length of the jets. The resulting prediction for the ratio of east jet to west jet 
brightness (for pieces of the jet of equal ages) as a function of the age of the components is shown in 
Figure O along with the observed intensity ratios. Second, the model Doppler factors D along the 
jet were used to predict the flux ratios. For a continuous jet the expected Doppler boost is D 2+a , 
while for discrete componen ts it is D 3+a , where a ~ 0.7 is the jet spectral index (for F cx u~ a ; 



Hjellming Johnston! ([1988)). The flux ratios predicted for both a continuous jet and for discrete 
components are also shown in Figured! We see that both projection effects and Doppler boosting 
vary slowly down the jet, and that the two are anti-correlelated. Neither projection effects nor 
Doppler boosting alone can account for the measurements; for example, the former predict that E2 
should be fainter than W2 and E3 fainter than W3, while the opposites are true. The combination 
of projection effects and Doppler boosting produces surprisingly good agreement with the observed 
ratios, especially for the model of a continuous jet. We regard this as evidence that the two jets 
are intrinsically very similar, and that both projection effects and Doppler boosts play roles in 
determining the observed intensity of the jets. More specifically, it confirms that the velocity of 
the jet material remains relativistic far down the jets from the optically-emitting region, consistent 
with what is required for the kinematic model to fit the larger-scale radio structure. Since (3 is 
relatively small in SS433, this analysis is not very sensitive to its exact value. 

Turning this analysis around, we can investigate the run of intrinsic brightness down the jet by 
removing both projection effects and Doppler boosting as described above. Observed flux densities 
and intrinsic brightness (in arbitrary units) are displayed against component ages in Figure [TUJ 
Consistent with the results above, we find that the two jets are intrinsically very similar, and 
that their brightness falls off with age (or distance from the core), roughly as a power law with 
/ oc i -1 oc r , or exponentially with characteristic time scale of about 75 days. 



In the conical sheath jet model of iHjellming &; Johnston! (|1988l ). jet intensity falls off with 



distance down the jet as / oc r n with n = (7p — l)/(6 + 65), where p is the electron energy 
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distribution exponent, and 5 = and 5 = 1 represent the freely-expanding and slowed-expansion 
cases, respectively. For a = 0.7, p = 2.4, and our mean index of n ~ 1 corresponds to 5 ~ 1.6, 
suggesting that the jets are significantly laterally confined as they propagate. 



3.2. Polarization 

The images in Figures [THSl show that the jets contain significant polarized fluxes at all frequen- 
cies. The fractional polarizations of the six jet components range from a few tenths of a percent at 
the three lower frequencies, up to m ~ 17% for E3 at 14.94 GHz. Along the jet but between the 
components, the fractional polarization is higher still, reaching m ~ 20% or more (see Q3. 41 below). 

The off-jet material is highly polarized; at the highest frequency, the fractional polarization 
is as great as m ~ 50 % (see Figure [5]) . This polarization cannot be due to electron scattering 
of photons from the core because the observed electric vectors are not circumferential. We are 
left to interpret the large m as due to synchrotron radiation where there is significant ordering of 
the magnetic field. This strong polarized signal may confuse the determination of magnetic field 
orientations in the jet itself (see §3.41 below). 

The small fractional polarization seen in the jets at the low frequencies is probably due either to 
side-side (beam) depolarization, or front-back depolarization consistent with large internal Faraday 
rotation. The east jet is significantly more highly polarized than the west jet at the low frequencies; 
this may be because the east jet is usually closer to us than is the west jet, and therefore suffers 
less depolarization near the source. 



3.3. Faraday Rotation 

Because SS433 has baryonic jets and is immersed in a stellar wind, we expect that Faraday 
rotation could be important. We examined this first by plotting the electric vector position angles 
X as a function of wavelength for each component at the three lower frequencies. Comparison of the 
middle panels of Figures HH3] shows a clear electric vector position angle rotation with frequency. In 
the case of a foreground Faraday-rotating screen, the observed electric vector position angles, the 
wavelengths, and the rotation measure RM are related by x = Xo + RM A 2 , where xo is the under- 
lying (unrotated) electric vector orientation. The results are consiste nt with an ave r age fo reground 



rotation measure of ~ 110 rad m 2 , in agreement with the results of IStirling et al.1 (12004! ) 



To examine the Faraday rotation distribution on smaller scales we derived rotation measures at 
each component using only the 4.80 GHz image and the convolved 14.94 GHz image. At frequencies 
of 4.80 and 14.94 GHz, a rotation measure of RM = 100 rad m~ 2 corresponds to rotations A^ of 
22° and 2.3°, respectively, and the 180° ambiguity in x corresponds to ARM = 898 rad m -2 . With 
only two frequencies, were are unable to determine unambiguous rotation measures; nonetheless, 
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we can estimate the size of the Faraday effect at each component by fitting the two frequencies to 
a A 2 -law, and choosing the position angle ambiguity that leads to the smallest RM. It is clear that 
the magnitude of the rotation measure at W3 is much larger than at the other components, and for 
it we determined the two smallest possible rotation measures magnitudes. This yields the rotation 
measures and intrinsic electric vector position angles in Table O 

A plot of Faraday rotation versus radial distance is shown in Figure [TT1 It is consistent with 
power law fall-off RM oc r~ n , with indices of n ~ 0.8 for the east jet and n ~ 2 for the west. 
Multi-frequency multi-epoch VLA observatio ns currently being r educed should show if and how 



this changes as the jets precess. We note that IStirling et al.l (|2004r ) also detected rotation measure 



gradients across SS433, and that their values are somewhat different than ours. The increase in 
rotation measure towards the core probably accounts for its negligible polarization through Faraday 
depolarization. 



3.4. Magnetic Fields 

Interpreting the radio emission as synchrotron radiation, the magnetic field vectors projected 
on the sky are perpendicular to the Faraday-rotation-corrected electric vectors xo- The distribution 
of inferred magnetic field vectors derived in this way is shown in Figure [12j Along the jets from E3 
to El and from W3 to Wl, the magnetic field directions curve gently with the jets, and lie roughly 
tangent to the ridge-lines. They are not align ed with the d i rectio n of ballistic motion of the jet 



material. This is consistent with the results of IStirling et al.l (|2004l ). The regions near El and Wl 



are harder to interepret because the jet trajectory is strongly curved and the polarized signal-to- 
noise ratio is rapidly decreasing. Note also that in the region E3-C-W3 the inferred magnetic field 
is roughly orthogonal to the jet locus. However, the fractional polarization there is small, and we 
do not regard this result as secure. 

The presence of highly-polar ized off-jet em i ssion is a surprise (high polarization of off-jet ma- 



terial is hinted at in figure 4 of IStirling et al.l ([2004)). In the regions north of E3 and south of 



W3 the magnetic field is remarkably unidirectional, lying at a position angle x — 45°. These high 
fractional polarizations show that the magnetic fields in the ambient medium are ordered on scales 
of thousands of AU. One possibility is that interactions between the jets and the ambient medium 
result in large scale motions of the surrounding gas that order an initially tangled field by shearing 
it. We do not know the three-dimensional distribution of this material, but it does appear to be 
trailing the jets, and thus probably lies inside the precession cone or on its boundary. 

The strongly polarized off-jet emission confuses the determination of the magnetic field orien- 
tations in the jet helices themselves. The jet ridge-lines are loci of minimum fractional polarization, 
with m higher in the adjacent off-jet regions. This suggests that the jet polarization is contami- 
nated by radiation from foreground and/or background ambient medium. We suggest that the jet 
magnetic fields are everywhere parallel to the observed helix. Between E3 and E2 and between 
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W3 and W2 the jet emission and the off-jet emission have the same magnetic field orientation. If 
the fractional polarization of the jet itself is smaller than that of the surrounding medium, a local 
minimum of m will be created along these parts of the jets. Between C and E3 and between C 
and W3, the jet and off-jet polarizations are orthogonal, with the off-jet radiation winning out, 
and the apparent polarization of the jet is lower than that of the adjacent emission. The polarized 
signals in the regions between E2 and El and between W2 and Wl are weaker, but they appear 
be consistent with magnetic fields aligned with the helix. 

At 14.94 GHz we can reliably determine the fractional polarization of the inter-component 
regions in the east jet, between El and E2 and between E2 and E3. From Figure [5] it is clear that 
m is larger in the these regions than in the adjacent components. In the west jet the same pattern 
occurs at lower signal-to-noise ratio. This suggests that the jets contain some tangled field that is 
partially ordered into a transverse orientation by compressions that cre ate the components . This 
process is apparently present in knots in parsec-scale jets of some AGN (IHughes et al.lll985l ). One 
extragalactic jet in wh ich the order in a predominantly longi tudinal field is apparently modified in 



this manner is 3C345 (Brown et al 



199 



Wardle et al.lll994l ) . If the components in SS433 are the 



loci of moving shocks, then some of the analysis in £13.11 might have to be modified; however, we 
know that the bulk velocity in the jets cannot vary significantly as they flow outward or else the 
kinematic model would fail to fit the structure beyond about 1". Only measurement of the proper 
motions of jet components can determine if they are moving at the jet fluid speed0 

Magnetic field lines that are tangent to the jet ridgeline and the strong leading-edge polarization 
suggest that the field orientations are determined at least in part by the interaction of the jet with 
the surrounding medium. For example, this is consistent with a model in which jet material plows 
into an ambient m edium conta ining tangled field, and the resulting compression orders the field 



in two dimensions (jLaing 



1981 



L|). Because of the precession of the jets, this process would occur 
predominantly on the "outsides" of the helices, and produce the sharp, highly-polarized leading 
edges that we see. Alternatively, shear created at jet-ambient medium interface could stretch a 
tangled field into a longitudinal one. 



4. CONCLUSIONS 

We have made polarization-sensitive VLA A-array images of SS4333 at five frequencies from 
1.4 to 15 GHz, using archival data over three days from March 1998. The principle results in this 
paper are: 

1. The structure of the SS433 jets at 4.80 and 14.94 GHz is consistent with the kinematic model 



1 Unfortunately, the ongoing upgrade of the VLA is resulting in EVLA antennas without 15 GHz receiver systems, 
a condition that will be remedied eventually. 
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and with a distance of 5.5 kpc if the velocities of ejection vary by about ±10% around the 
optically determined value. 

2. The 14.94 GHz image of the jets reveals a continuous jet, on top of which there are four pairs 
of brighter components. Comparison with the kinematic model shows that the two members 
of each pair were emitted simultaneously and with equal velocities. Pairs emerge at intervals 
of about 35 days, but the process is not strictly periodic. 

3. When the observed brightness ratios of the pairs of components are corrected for projection 
effects and Doppler boosting, we find that the two jets are intrinsically very similar. 

4. The intrinsic intensity of the jet components falls of with their age (or distance from the core) 
roughly as a power law with index —1, or exponentially with a time scale of about 75 days. 

5. The jets are linearly polarized with fractional polarization reaching m ~ 20%, consistent with 
synchrotron radiation from relativistic electrons. The core is essentially unpolarized, probably 
due to Faraday depolarization. 

6. The Faraday rotation is consistent with a foreground screen with RM ~ 110 rad m -2 , plus a 
larger contribution of several hundred rad m -2 close to the core that is asymmetrical across 
the source. 

7. The jets are surrounded by significant off-jet emission that is primarily trailing the jet ridge- 
line. This emission is highly linearly polarized, ranging as high as m ~ 50% or more, indicating 
the presence of large-scale magnetic fields and relativistic electrons in the ambient medium. 
The jet ridgeline is a locus of minimum fractional polarization, suggesting contamination by 
foreground and/or background off-jet emission, and confusing the determination of the field 
orientation in the helices themselves. Taking this into account, the polarization images are 
consistent with a magnetic field that is everywhere parallel to the jet locus. 

8. The fractional polarization of the jets is greater between the components than in the compo- 
nents, suggesting that the jet contains some tangled magnetic field that is partially ordered 
by compressions at the components, as seen in extragalactic jets. 

9. The fact that the leading edges of the jets are sharper in profile than the trailing edges and 
the longitudinal orientation of the magnetic field show that significant interactions between 
the jets and the surrounding medium play an essential role in determining the morphology 
and magnetic field of SS433 jets on the scales seen in VLA images. 



Further polarization-sensitive observations covering a range of precession phases and at multiple 
frequencies that are currently being analyzed should shed more light on the role of jet-ambient 
medium interactions in SS 433. This may lead to new perspectives on similar phenomena taking 
place in the parsec-scale jets of AGN. 
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Table 1. Ephemeris for SS433. 



Property 


Symbol a 


Value 


Reference 13 


Jet Speed 


13 = v/c 


0.26 ± 10% 


1 


Precession Period 


P 


162.375 d 


2 


Reference Epoch (TJD) 


tref 


3563.23 


2 


Precession Cone Opening Angle 


ij, 


20.92° 


2 


Precession Axis Inclination 


i 


78.05° 


2 


Precession Axis Position Angle 


X + vr/2 


98.2° 


3 


Sense of Precession 


Srot 


-1 


4 


Distance 


d 


5.5 kpc 


5 



a Following iHjellming Johnston! (|1981bl ) . 
b References-(lj See tex t (2) lEikenberrv et al I d200lh . 



Mioduszewski 



( 200m. ( 4) iHiellming fc Johnston 



Blundell fc Bowleri feooj ) ; lLockman. Bhmdell. fc Gossl hooJo . 



981b), 



(3) 
(5) 
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Table 2. Components in the Jets of SS 433, Epoch 1998 March 5-7. 



z^(GHz) J(mJy) p(mJy) m{%) \ 



El: (+1.008, +0.153)", (1.020, 81°) fe , r = 121 c 



4.800 


16.5 


1.83 


11.1 


+9' 


14.94 


2.16 


0.140 


6.5 


-22' 


14.94 d 


4.68 


0.225 


4.8 


-8' 



E2: (+0.745, +0.061) a , (0.747, 85°) b , r = 80 c 



4.800 


34.1 


4.93 


14.8 


-r 


14.94 


4.24 


0.245 


5.8 


-7 


14.94 d 


10.4 


1.08 


10.4 


-30' 



E3: (+0.450, -0.152) a , (0.475, 109°) 6 , r = 53 c 

4.800 59.9 6.39 10.7 -7° 
14.94 8.31 1.43 17.2 -37° 

14.94 d 20.1 2.55 12.7 -42° 



C: (0,0) a , (0,0) 



4.800 


398. 


0.840 


0.2 


-55' 


14.94 


196. 


0.349 


0.2 


-15' 


14.94 d 


214. 


0.518 


0.2 


-47' 



W3: (-0.331, +0.135) a , (0.357, -68°) fe , r = 48 c 



4.800 


78.0 


1.31 


1.7 


+70° 


14.94 


4.72 


0.431 


9.1 


-46° 


14.94 d 


26.2 


0.745 


2.8 


-49° 



W2: (-0.610, -0.073) a , (0.614, -97°) fe , r = 83 c 
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Table 2 — Continued 



i/(GHz) 


J(mJy) p(mJy) m(%) 


X 


4.800 


21.6 0.294 1.4 


+31° 


14.94 


2.14 0.217 10.1 


-4° 


U.M d 


5.88 0.390 6.6 


-16° 


Wl: (- 


0.934, -0.151) a , (0.946, -99°) 6 , 


r = 120 c 


4.800 


10.8 1.08 10.0 


+20° 


14.94 


1.19 0.148 12.4 


-2° 


14.94 d 


2.41 0.128 5.3 


+ 3° 



a Position relative to the core, as (a, 5) in arcsec- 
onds. 

b Distance from the core in arcseconds and struc- 
tural position angle. 

Mo del component age in days. 

d Values measured at the resolution of the 
4.800 GHz images. 
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Table 3. Faraday Rotation in SS 433, Epoch 1998 March 5-7. 



Component 


RM (rad m~ 2 ) 


Xo 


9 B a 


El 


+88 


-11° 


+79° 


E2 


+144 


-33° 


+57° 


E3 


+173 


-46° 


+44° 


W3 b 


-302 


-42° 


+48° 


W3 6 


+596 


-63° 


+27° 


W2 


+235 


-22° 


+68° 


Wl 


+84 


+1° 


-89° 



a Position angle of the underlying magnetic 
field. 

b Two possible rotation measures. 
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RcJaLivc RighL Ascension (arcscc) Relative RighL Ascension (iarcscc) 

Fig. 1. — VLA images of SS433 at 1.415 GHz, made from data taken over 1998 March 5-7. 
(Top) Contours of total intensity, starting at 0.75 mJy beam -1 , with steps of factors of v2, and 
peak of 811 mJy beam -1 . (Middle) Linear polarization, showing contours of polarized intensity, 
starting at 0.75 mJy beam -1 , with steps of factors of \/2, and peak of 10 mJy beam -1 , with 
ticks indicating electric vector position angles, plus a single total intensity contour. (Bottom) 
Fractional polarization on the scale to the right, overlying every-other total intensity contour. The 
dots indicate the locations of seven bright total intensity components, as determined at 14.94 GHz 
(Figure [5]). The kinematic model described in the text is shown as material ejected at five day 
intervals, with "+" and "o" indicating oncoming and receding parts, respectively. The restoring 
beam with FWHM of 1.38" x 1.38" is shown as a cross to the lower-left. (Both monochrome and 
color figures are included for monochrome printers.) 
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'10 5 -5 -10 '10 5 -5 -10 



RcJaLivc RighL Ascension (arcsec) RclaLivc RighL Ascension (arcsec) 

Fig. 2. — VLA images of SS433 at 1.465 GHz, made from data taken over 1998 March 6-7. Same 
parts as Figure[Q The contours of total intensity start at 0.75 mJy beam -1 , have steps of factors of 
\/2, and a peak of 732 mJy beam -1 . The contours of polarized intensity start at 0.75 mJy beam -1 , 
have steps of factors of v2, and a peak of 5.9 mJy beam -1 . The restoring beam with FWHM of 
1.10" x 1.10" is shown as a cross to the lower- left. (Both monochrome and color figures are included 
for monochrome printers.) 
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'10 5 -5 -10 '10 5 -5 -10 



RcJaLivc RighL Ascension (arcsec) Relative RighL Ascension (arcscc) 

Fig. 3. — VLA images of SS433 at 1.665 GHz, made from data taken 1998 March 5. Same parts 
as Figure [TJ The contours of total intensity start at 0.75 mJy beam -1 , have steps of factors of \^2, 
and a peak of 775 mJy beam -1 . The contours of polarized intensity start at 0.75 mJy beam -1 , 
have steps of factors of y/2, and a peak of 5.5 mJy beam -1 . The restoring beam with FWHM of 
1.09" x 1.09" is shown as a cross to the lower-left. (Both monochrome and color figures are included 
for monochrome printers.) 
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2 -2 2 



RelaLivc Right Ascension (arcaoc) liclaLivc Right Ascension (arcacc) 

Fig. 4. — VLA images of SS433 at 4.80 GHz, made from data taken 1998 March 5-7. Same parts 
as Figure [H The contours of total intensity start at 0.20 mJy beam -1 , have steps of factors of y/2, 
and a peak of 398 mJy beam -1 . The contours of polarized intensity start at 0.20 mJy beam -1 , 
have steps of factors of \/2, and a peak of 8.1 mJy beam -1 . The restoring beam with FWHM of 
0.39" x 0.39" is shown as a cross to the lower-left. (Both monochrome and color figures are included 
for monochrome printers.) 
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RelaLivc Right Ascension (arcaoc) liclaLivc Right Ascension (arcacc) 

Fig. 5. — VLA images of SS433 at 14.94 GHz, made from data taken over 1998 March 5-7. Same 
parts as Figured! The contours of total intensity start at 0.15 mJy beam -1 , have steps of factors of 
\/2, and a peak of 196 mJy beam -1 . The contours of polarized intensity start at 0.15 mJy beam -1 , 
have steps of factors of v2, and a peak of 1.5 mJy beam -1 . The restoring beam with FWHM 
of 0.125" x 0.125" is shown as a cross to the lower-left. (Both monochrome and color figures are 
included for monochrome printers.) 
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RelaLive Right Ascension (arcsec) Relative Right Ascension (arcsec) 

Fig. 6. — VLA images of SS433 at (a) 4.80 GHz and (b) 14.94 GHz, compared to three kinematic 
models. In each image we plot models with (top) the optical velocity of (3 = 0.2647 and distance 
of d = 5.0 kpc, (middle) jet with (5 = 0.2647 and d = 5.5 kpc, and (bottom) jet with (3 = 0.29 and 
d = 5.5 kpc. 
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1.5 1 0.5 0.5 -1 -1.5 



Relative Right Ascension (arcscc) 

Fig. 7— VLA image of SS 433 at 14.94 GHz, made from data taken over 1998 March 5-7, with 
the core removed. The contours of total intensity start at 0.15 mJy beam -1 , have steps of factors 
of \/2, and a peak of 15 mJy beam -1 . The dots indicate the locations of six jet components, as 
determined from Figure 02 the tentative components E4 and W4 found here are also labeled. The 
kinematic model described in the text is shown as material ejected at five day intervals, with "+" 
and "o" indicating oncoming and receding material, respectively. The restoring beam with FWHM 
of 0.156" x 0.129" is shown as a cross to the lower-left. 
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East Component Age (days) 



Fig. 8. — The model ages of the four pairs of jet components in SS 433, plotted as age of west compo- 
nent versus age of corresponding east component. The innermost point is not used in the statistical 
analysis. The sloping line represents the hypothesis that the pairs were emitted simultaneously. 
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Fig. 9. — Predicted versus measured flux ratios in SS 433 at 14.94 GHz for four pairs of components 
in the jet of SS433, plotted as a function of the average age of each pair of components. Model 
predictions are based on d = 5.5 kpc, v = 0.29c, and a = 0.7, illustrating (a) projection effects 
(dotted line), (b) Doppler boosting for a continuous jet (lower dashed line), Doppler boosting for 
discrete components (upper upper dashed line), and (c) these two effects combined (lower and 
upper solid lines, respectively). 
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SS 433 I99K March 5-7 
14.94 GHz v = 0.29c, (1 = 0.7 
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Component Age (days) 



Componcnl Age (tlays) 



Fig. 10. — (a) Observed and (b) intrinsic jet fluxes in SS433 at 14.94 GHz for four pairs of 
components in the jet of SS 433, plotted as a function of the age of the components. In both parts 
(a) and (b), components in the east jet are shown as open symbols, those in the west jet as filled 
symbols. In (b), models for a continuous jet are shown as circles, those for isolated components as 
squares. 
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Fig. 11. — Rotation measures in SS433 as a function of radial distance from the core, for the six 
bright jet components described in the text. Shown is the rotation measure determined from the 
4.80 and convolved 14.94 GHz images. The power law fits are described in the text. 
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Relative Right Ascension (arcsec) 



Fig. 12. — Magnetic field orientations in SS433. The contours are polarization intensity, with tick 
marks of Faraday-rotation-corrected magnetic field vectors, and a single total intensity contour at 
0.15 mJy beam -1 . The polarization contours start at 0.15 mJy beam -1 , have steps of factors of 
\/2, and a peak of 1.5 mJy beam -1 . The kinematic model described in the text is shown as material 
ejected at five day intervals, with "+" and "o" indicating oncoming and receding parts, respectively. 
The restoring beam with FWHM of 0.125" x 0.125" is shown as a cross to the lower-left. 



